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ABSTRACT 
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was designed, built, and tested. A calculator program was 
written which enabled a research group studying laser 
propagation in the marine boundary layer to conduct rapid, 
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I. INTRODUCTION 


The potential utility of laser technology in the areas 
of line-of-sight communication, target tracking, and wea- 
ponry makes an understanding of the propagation of laser 
beams near the air-ocean interface imperative. The effects 
of a turbulent atmosphere on laser transmission were theor- 
etically modelled by Tatarski [Ref. 1] et. al. and, subse- 
quently, extensive field research in the terrestrial boundary 
layer has substantiated the bulk of the theory. Research 
data in the marine boundary layer has been scarce and 
inconclusive. 

Atmospheric effects on laser propagation may be divided 
between laser-induced phenomena and those caused by meteoro- 
logical characteristics independent of the laser's presence. 
The first classification applies only to lasers of sufficient 
power to heat the atmosphere (thermal blooming, upwind 
drift). The second is a universal grouping affecting the 
whole electromagnetic spectrum regardless of power. The 
atmospheric effects include beam refraction, beam spread 
(or degradation of resolution), absorption, and scattering. 
Refraction may be further categorized according to the 
dimensions of the atmospheric whirls. Large-scale vortices 
cause path deviation (beam wander) while small-scale vortices 
cause time-intensity variations (scintillation) and degrada- 


tion. oO: resolution. 





The Naval Postgraduate School, Monterey, has been involved 
over the last four years in a continuing effort to establish 
the data base necessary to confirm or modify existing theory 
as applied to the air-ocean interface. Two research groups 
have operated simultaneously but independently; the optical 
research group and the atmospheric turbulence group have 


attempted to correlate determinations of C,, (the refractive 


N 
index structure constant) through the respective analysis 

- light intensity distributions and meteorological data. 

The optical research group has concurrently investigated 
scintillation, extinction (scattering and absorption) and 

the atmospheric modulation transfer function (a measure 

of beam spread and wander). Theses on the effects of aperture 
averaging, determination of an atmospheric extinction con- 
stant, and determination of an atmospheric MTF for a black- 
body source were completed over the same period as the 
agener s {[Refs. 10, li, 12]. 

A major problem for the optical group was the inordinate 
amount of time necessary to process scintillation and extinc- 
tion data. The methods employed depended heavily on operator 
analysis of raw data or extensive manual input to computers. 
While the processing theory being applied was appropriate, 
the time penalty and the inherent high probability of error 
had to be minimized. The author concentrated his efforts 
in this area employing only the basic results of laser 


propagation theory. 





If. PROJECT BACKGROUND 


Various members of the Naval Postgraduate School's 
faculty have been conducting research on atmospheric pro- 
pagation phenomena since 1972. The first experiments inves- 
tigated the atmospheric MTF and scintillation in the basement 
corridor of Spanagel Hall. In the scintillation experiment 
a laser beam reflected from a flat mirror placed at the 
Opposite end of the 137 meter passageway and entered a pin 
diode detector located near the laser source. The detector 
Output, proportional to the incident beam intensity, was 
passed through a log-amplifier and recorded on magnetic tape. 
The data was later periodically sampled, digitized, and 
accumulated by a pulse-height analyzer; the resultant dis- 
tribution was Gaussian, as theory had predicted. Refractive 
index structure constants were calculated from photographs 
of the distributions by measuring the curves’ full width at 
half maximum height. The MTF measurement used a six inch 
parabolic mirror mounted as a Newtonian telescope. The 
image from a laser at the other end of the passageway was 
scanned across a slit using a cam-driven rocking diagonal 
mirror. The output of a photocell was fed to a boxcar inte- 
grator to generate time-averaged data, which was then recorded. 
Turbulence in the hall was somewhat controllable and detec- 


table reduction of MTF with turbulence was observed. 
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In late 1972 a group of interested researchers discussed 
possible environmental experiments. The first such experi- 
ment was completed 10 January 1973 by faculty, staff, and 
students and involved shining a low-power (mW) red laser 
(6328 A) across a 4.3 kilometer range at the southern end 
of Monterey Bay. The research vessel Acania was anchored 
at the approximate midpoint of the beam path to take atmos- 
pheric turbulence measurements with onboard sensors. Ina 
second experiment at San Nicholas Island (3 March 1973) low 
power laser beams were transmitted from the northwest tip 
of the island to the Acania anchored offshore. The qualified 
success of these experiments completed with equipment already 
owned by the school was encouraging enough to prompt the 
submission of a research proposal to the Naval Ordnance 
Laboratory (White Oak, Md.) for fiscal year 1974. 

Experimental apparatus has changed significantly since 
the Navy lab started sponsoring the research. A Cassegrainian 
telescope system with an eighteen inch parabolic mirror and 
an elastically mounted scanning mirror (driven electromag- 
netically by a shaped waveform) was used in the most recent 
MTF measurements. An enclosed trailer was built to transport 
and protect the telescope and associated equipment; aluminum 
support legs that extended to rest directly on the ground 
negated the previously persistent problem of telescope motion 
due to operator movement. More recently a bus has been used 


which holds most of the on-shore optical and meteorological 


ae 
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equipment and employs the same independent telescope suspen- 
sion idea. Ship motion originally made it difficult to keep 
even a diverged laser beam on the telescope long enough to 
collect maningful MTF data. A gyro-stabilized tracking system 
was designed to go aboard the Acania. At first it was manned 
by two operators who sight-aimed a gyro-mounted laser with 
hand-operated servo joy sticks; one controlled horizontal 

and the other the vertical train. Little useable data was 
obtained with this system because of problems with human 
coordination and fatigue. The solution was a quadrant detec- 
tor servo system which will lock on a laser beam once the 
operator has aligned it with the target. Shipboard scin- 
tillation detectors are also mounted on the gyro frame and 
thereby maintain a constant relative position in the incident 
beam path. The gyro has been housed in an aluminum shed, 

when aboard the Acania, to protect both equipment and operator 
from the elements. 

Scintillation data was best taken aboard the Acania; a 
detector with a wide angle of acceptance could continuously 
sample a diverged beam from ashore more readily than the 
reverse because of the effect ship motion would have on the 
direction of propagation of a beam originating on board. 

The single-detector system originally employed required 
excessive attention during data runs because light from 
environmental sources refracts into the beam path and can't 


be distinguished from that of the monitored laser. Since all 
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photocell data passed through a log-amplifier a skewed 
distribution would result. In an attempt to negate this 
effect the beam had to be interrupted ten seconds out of 
every minute to zero-out detector voltage generated by back- 
ground radiation. A dual-detector system was built with the 
idea that the output of one pin diode detecting only background 
could be subtracted continuously from that of another unit 
aimed at the laser beam. 

Since the detectors (supposedly identical) produced 
different proportional outputs extensive voltage biasing 
was necessary. An even greater problem was that ship motion 
continuously changed the position of the detector relative 
to the beam center. Subsequently another single-detector 
system was designed, this one employing a chopped laser source 
and a demodulator. The chopping was accomplished at ane source 
By reflecting the laser light first off of a cylindrical 
mirror which fanned the beam horizontally and then off of 
a flat mirror vibrating at approximately two kilohertz which 
made the beam scan vertically. The demodulator unit took 
the detector output and electronically subtracted the minima 
(background) from the maxima (laser & background) and sent 
the difference to the log-converter. Detector sensitivity 
was recently increased by replacing the pin diode with an 
avalanche detector. An FET follower converts the high- 
impedance output of the detector to a low impedance which 


can be transmitted a greater distance through coaxial cable. 


EG 





This has allowed all data processing gear on the Acania to 
be positioned below-decks in the dry lab. 

Processing of scintillation data was initially best 
accomplished through physical measurement of photographs 
of resultant distributions displayed on the PIP-400 pulse 
height analyzer. Attempts to record and subsequently analyze 
data with integral methods on a computer suffered from poor 
Pegtei zing. Digital Gutput could be obtained directly from 
the PIP-400 by a fast printer or from a Nuclear Data 128- 
channel analyzer by teletype. An attempt was made to fita 
Gaussian curve to the data using least-squares method. This 
analysis proved highly successful when extraneous voltages 
did not skew the experimental distribution but was painfully 
slow because data had to be manually input to the computer 
point by point. The most recent improvement has been the 
direct interfacing of the PIP-400 with a compact computer 
programmed to process scintillation and extinction data on- 
scene within six minutes of sampling. This was completed 


by the author and is the subject of this thesis. 
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A. SCINTILLATION 

The theory of optical propagation in a turbulent atmos- 
phere has been extensively documented and RAbee portions 
employed by the Naval. Postgraduate School optical propaga- 
tion group have not changed significantly since the project 
began. A thesis submitted by B. C. Haagensen [Ref. 3] in 
June, 1973 contained an excellent outline of the pertinent 
theoretical derivations and an extensive bibliography. The 
author does not intend to list or verify familiar theoretical 
developments but rather to present the significant end 
products vital to data processing. 

Theory states that atmospheric turbulence is a random, 
non-linear, three-dimensional phenomenon. A turbule is 
ideally assumed to be an isotropic, homogeneous pocket or 
Swirl of air with a particular index of refraction. The 
inertial subrange delineates a region of turbule size in 
which refraction is theoretically described; it is bounded 
on one side by turbules so small (a few millimeters) that 
viscous dissipation of energy causes their immediate collapse 
and on the other by eddies so large (order of meters) that 
they are no longer homogeneous and energy is added to the 
turbulence spectrum. In the inertial subrange large turbules 
form and transfer energy to successively smaller turbules 


with minimal losses through viscous dissipation. 


dee, 





Tatarski [Ref. 1] described the three-dimensional power 


spectral density of turbulence, o (K), with the equation: 


where: 


c 


N 


“K*, 
Bi 


m 


= 0.033 ¢.” K 77/9 exp | 


refractive index structure constant 


au 


Spatial wave number = ? 


coherence length of a turbule 


oat Lo = lower limit of inertial subrange 


obviously parameterizes the total amount of energy 


in the turbulence. The random distribution of turbulence 


cannot be forecast and must be described statistically. A 


beam is modulated in a multiplicative manner when it passes 


through successive turbules of differing index of refraction. 


It follows that the natural log of the beam intensity is 


altered additively and should follow a normal probability 


distribution as a result of the central limit theorem. Cy 


1s directly related to the structure function of the refrac- 


tive index (the mean square difference of the index of 


refraction at two points a distance d apart) by: 


= <[n(r) - n(r+d) ]°> = eve qz/3 
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After lengthy derivation and a number of assumptions [Ref. 
oy egw 10), Cy can be seen to be directly dependent on the 


variance of the log intensity distribution by: 
2 Silly Vig. ea Ai ee 


N Cert 


where 


<finthe— haeeses = variance 


Q 
fo 
IH 

lI 


N 
lI 


optical propagation path length 


Cy can also be calculated by a statistical determination 


of the temperature structure constant, C Since tempera- 


T é 


ture is the dominant factor ina meteorological formulation 


of the index of refraction: 


_ P 8 
Cy = (79.0 72 cn) Cm 
where 
P = mean pressure in millibars 
Ih = mean temperature in K 
2 
C Me -Gueee = Ti rtd) } > 
at 2/8 


d 


2 





A major objective of the project is to study the 
correlation between Cy values obtained optically and 


meteorologically. 


B. GAUSSIAN MODELLING 

According to theory, periodic sampling of a laser beam 
propagating in the inertial subrange should approximate a 
lognormal distribution; if the natural log of the detector 
voltage (proportional to beam intensity) defines abscissa 
values and the ordinate marks the number of samples returning 


a specific voltage the resultant plot should be of the form: 


(x ~ ai 
20 


y = C exp [- 


where: 


C is a constant 


x is natural log of detector voltage 


» | 


is the average natural log of detector voltage 


o* is the variance i) 


Taking the natural log of each side: 


2 = —2 
iia. =. — -— Ts _ + (tree 7 = ax* + Dea ae 
oj 


Zo 20 


Ze 





where: 


es 
Ee 
_ 
a ae 
0} 
<2 
c=tlincere—_~_ 
20° 


Experimentally the natural log of detector voltage axis 
is divided by the PIP-400 into two hundred channels of equal 
width and the y axis is the number of counts (or times) the 
equivalent detector voltage (or beam intensity) was detected 
over the sampling period. The number of counts in each chan- 
nel after a data run can be assumed to be Poisson-distributed; 
for one measurement the 68% confidence interval is y + Vy. 
Since the natural log of y Ss Helen to obtain the quadratic 


form, the uncertainty in ln(y) becomes: 


Siemon] - = 


[+ 


So, for the log distribution, the channels with the largest 
numbers of counts have the narrowest confidence intervals 
(smallest uncertainty) and should be weighted more heavily. 
The effect of narrowing confidence intervals with increasing 
counts is easily illustrated by plotting them on log paper. 

In the original program it had been decided to weight the data 


inversely with the uncertainty, i.e. multiply by Vy. 


ce 





THE LOGNORMAL 


LN X ? 








The ln(y) distribution resulting experimentally can be 
approximated with a polynomial by means of least-squares 


fitting. The method seeks to minimize 


eels, “ 5 
i w, [in(y,) - ax, 9 - Dx, - cl 
where: 
x, = channel number 1 
en counts in channel i 
a a weighting of data in channel i, VY, 


Setting the partial derivatives with respect to a, b, and 
c equal to zero and cancelling common terms gives three 


equations with three unknowns: 


2 _ 4 3 2 
Iw, X, iny,; = alw,x, + biw.X; 3 Clw, xX, 


a 3 2 
iw,xX,iny, = alw. xX. + biw, xX, + CLw. xX, 


= 2 
iw, iny, = alw,X, + biw. x, sr CLw, 


Applying Cramer's Rule [Ref. 6] to the system determines 


the coefficients: 
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det (A, ) 


a * det (A) 
det (A,) 
= 
det (A) 
det (A) 
=— det (A) 
where 
det(A) = Sw. xX. “Ew x. “{rw = EW eee 
suaeak eo 1 11 
- TW Xs "tow, x ow, 5G. ,7mW, age Sed 
wheel aL 
+ IW Xs * tow, ae lw, mee gg OH oe 
= 3 = 
det (A,) = Iw. x, {lw, x, Zw, iny, Zw, lw,x,lny,} 
+rw.xX 2 sus on Wee thy = =). W > X ier xX einy } 
cae a saa eats eae ie deers! ack 1 
Z 2 
+ Iw, x, {lw, Xs TW See ,iny; —Iw; Xs TWX, lay; } 
= 4 = 
det (A,) = Iw, xX, (lw, lw, x,lny, Iw. X,2w,iny, } 


+ Dw .x,~{2w, m. “aw, ,iny; ~ tw UW Ge polnyee 
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_ 4 2 
det (A,) LW x, {lw x, tw,iny, Uw, XW. x, 


Z 


3 : a 3 
tiw. x, {iw,xX,Zw, x, any ,-lw, x, tw. lny, } 


2 3 - 2 2 
+ Ew. x, tlw, x, lw, x, lny,-Zw, x, “Lw.x, lny, } 


Exactly the same forms for the coefficients result from 
applying Gauss-Jordan elimination [Ref. 6]. A total of 
eight distinct weighted-data summations form the basic ele- 
ments which combine to give the coefficients. Coefficient 
"a" and a channel number-to-log voltage conversion factor 
deliver the lognormal variance ee) which, with the optical 
wavelength and range, determines the refractive index struc- 
nT constant (C,). Sceeticientseta', "bl", and the conversion 


factor fix the average log voltage value (ln 


eae LATINCTION 

Extinction is the loss of time-averaged intensity with 
distance by a beam propagating through the atmosphere. No 
extinction measurements had been taken by the optical propa- 
gation group before September, 1975. Assuming an exponential 
dependence and spherical spreading, extinction can be 


represented by an equation of the form: 





— 0 = a-_ 
I= 5 spo [Sans Ry) J 
where: 


Za 





I = intensity at distance R 


Ty = intensity at distance R 


0 
Ro = some arbitrary standard distance 
R = some arbitrary distance 
uw = the extinction ccefficient 


An equivalent form of the equation is: 
I = Sir R * exp (uR )] exp(-uR) 
R- 0-0 0 


If R, is some small distance (on the order of a meter), the 


0 
exponential with Ro approaches one and the above form can 


be simplified EO: 
ae 
1 exp (-yR) 
R 
where 
Ca=) 1 -R (a measure of source intensity) 
Taking the natural log of both sides: 
inI = lnc - 21lnR - uR 


Solving LOr i: 
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Ln (—S5) 


= 2n¢ = 2inR = int _ __iR]_ 
R R 


The constant C must be determined at the time of the experi- 
ment and the time-averaged intensity and range must be 
measured. Since detectors generate voltages proportional to 
the incident intensities, voltage can be substituted for 
intensity with meaningful results. Phil Parrish investigated 
the extinction constant (C) and coefficient (u) in a con- 
current thesis [Ref. 10]. 

A statistical analysis was done to determine the extinc- 
tion coefficient's degree of sensitivity to seo. Ln ee 
Ree and L: 


nen.) inc - 2inR - lnI 


Ju _ 172 


= f (ou Z 2 Ol 2 Z 


du, 2 2 
sk)“ act + (g4)* ar® + (gE)* ar*} 


oe] 


Seecking Sensitivity to error with the true value of 1 


assumed to be l x La: 
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Obviously, errors in all quantities have more effect 
at shorter ranges. Errors in the extinction coefficient 
appear to depend linearly on independent errors in the 
extinction constant, time-average intensity, and range but 
decrease with increasing true range. Errors in range have 
a much greater effect than equivalent errors in either of 
the other quantities. The last two calculations in the 
previous table use maximum error values estimated to be 


appropriate with present equipment and procedure. 
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IV. INTERFACING 


iio wicrlnguon Normal distributions to scintillation 
data through a computer program using least-squares method 
had proven highly successful but also extremely slow. The 
PIP-400 pulse height analyzer used to accumulate data dis- 
tributions had the design capability to output the information 
through a 50-pin plug in parallel BCD (binary coded decimal) 
readout. The HP9810A calculator in which the processing took 
place had an optional parallel input interface card which 
accepted BCD input through a 50-pin plug at the same voltage 
levels as put out by the PIP-400. It was thought that a 
Simple interface matching the proper lines would obviate the 
pokey manual input of data. A number of hidden problems 
were found which made the project much more difficult than 
expected. 

The first problem was that both machines had been designed 
as masters. The analyzer sent a high-to-low (voltage) print 
command and waited for a low-to-high flag from the peripheral 
(which told it to step to the next channel). The calculator 
sent a high-to-low voltage command to an external source 
asking for data and waited for a low-to-high flag which told 
it the information was ready. The solution was to invert 
both sets of signals so that the data request from the calcu- 
lator appeared as a peripheral flag to the analyzer anda 


data print command from the PIP-400 appeared as an information 
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ready flag to the HP9810. If the analyzer function knob 

is set to “PRINT” before the calculator asks for data the 
first channel will be ignored, but since the first (and 
last) five channels of the distribution are zeroed by the 
program anyway (because of the high probability of extrane- 
ous counts in these channels) the deletion has negligible 
effect. Since the voltage signal levels for both machines 
matched, the modification was completed with a single hex 
inverter chip (Signetics 7404) using two of the six inverters 
available. A six volt supply lead from the analyzer and a 
100-ohm resistor filled the inverter's power requirement 
won) 25> volts). 

Matching the data lines proved more of a problem. Data 
from the analyzer in parallel readout consists of six data 
digits followed by three address digits. Digital information 
is relayed in BCD which requires four bitlines for each 
decimal digit. The PIP-400 puts out a low current, high 
impedance data signal while the HP9810 uses a pull-up resis- 
tor to bring an expected high current, low impedance signal 
mp to five volts for a "1" bit. The calculator reads any- 
thing over 2.8 volts as a "1" bit and anything less as a 
men Dit. With the command lines being inverted by the chip 
and the thirty-four information bit lines properly mated 
(inside a junction box with two 50-pin plugs) a calculator 
print-out of analyzer data showed all 9's; this meant that 


the pull-up resistor in the calculator was raising all 
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analyzer output signals to over 2.8 volts. This supposition 
was substantiated when the correct data was properly read 

out of the analyzer by a fast printer through the same output 
mud. 

The first solution attempted was to install a grounded 
quarter-watt resistor in parallel with each signal line to 
pull the voltage of a "0" bit below 2.8 volts. The address 
numbers were successfully printed using a 5600 ohm resistor 
on each lead but the data digits were decidedly non-uniform 
in their individual bit resistor requirements. Comparison 
of calculator printouts with fast printer readouts were used 
to determine the necessary resistance; lead requirements 
varied between 4700 and 8200 ohms for the calculator to print 
data digits accurately. Variable (0 - 10,000 ohm) poten- 
tiometers were used to find the range of resistance values 
between those which fixed a bit's output at either "0" or 
"1". The proper resistance span was smaller, in some cases, 
than that separating succesSive standard resistors and was 
found to vary with time (possibly temperature dependent). 
Upon testing line voltage levels it was found that while the 
BoGEctorns broward) O° bit £rom 3.0 to 2.7 volts it also 
scale an i = btemcdovierrom 5.0 to 2.9 volts (approximate 
values). Since the fixed-resistor interface was obviously 
too finicky it was abandoned. 

The next scheme attempted employed operational amplifiers 


(Signetics A747) constructed on silicon chips as unity-gain 
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voltage followers and succeeded in generating the high 
current, low impedance input the calculator was designed 

to use. Seventeen op-amps were required to handle the thirty- 
four bits of information which made up each data point. The 
chips were all connected in parallel with a 100 mA dual power 
Supply (run on 115 volt AC) which generated the + 15 volts 


direct current necessary. 


V = -15 Volts 


Operational amplifier 
socket wiring for use 
as a voltage follower 


input A 
Output A 
(A & B are BCD bit Lines) 
Output B 
input B 





V = +15 Volts 


With the op-amps and additional wiring the circuitry was 


too bulky to be contained in the original junction box; after 
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the system had proven itself in the field it was redesigned 
and built with a junction box large enough to hold all the 
components internally. 

The interfacing of analyzer and calculator cut the time 
to input each data point from seconds to milli-seconds; 
with a 200-channel analyzer (PIP-400) this meant a savings 
of five to ten minutes processing time on each run and made 
timely on-scene analysis of experimental distributions 


possible. 
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The remnants of the first interface uSing operational 
amplifiers. While successful it was too vulnerable in 
an operational environment. 





The interface in its final form. To the right are the op- 
amps and inverter chip, to the left is the power supply. 
Deemieelug om rignt goes to calculator, one on left to 
analyzer; outlet plug is for dual power supply. 
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V. COMPUTER PROGRAM DEVELOPMENT 


A. MODIFICATION OF THE ORIGINAL PROGRAM 

While mechanically interfacing the pulse height analyzer 
(PIP-400) and the computer (HP9810) saved a great amount of 
time previously used to input data, the computer program 
which processed the data was itself inefficient. An inor- 
dinate amount of time was spent labelling, giving operator 
instructions, and searching for subroutines. The program 
had been written in a straight-line manner with a minimal 
amount of branching for recurrent processes. The original 
program processed data for a least-squares fitting of a 
Gaussian and calculated the lognormal variance (SIGMA) and 
the refractive index structure constant (Cx) bute did not 
compute the average detector voltage or an extinction coef- 
ficient. The program lacked the versatility to accept data 
from available analyzers by manual input. 

The Hewlett-Packard 9810A calculator (with option 003 
installed) has a total of 2036 locations in memory available 
for program steps of which 1847 were used in the original 
program. Since the additional desired processing could not 
possibly be accomplished in the steps remaining, a major 
program overhaul was required. The first measure taken was 
to delete or shorten as many labels and operator instructions 
as possible; approximately 250 program steps were excised 


in this manner with no functional loss. More steps were 
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saved by recording implemental constants on a magnetic card 
and entering them directly into the proper memory locations 
immediately after the program was loaded. 

The program now had enough room for additional processing. 
In previous analysis [Ref. 3, et.al.) it was assumed that 
the natural log of the average detector voltage (beam inten- 
Sity) was approximately equal to the average natural log of 
detector voltage. The validity of this assumption, initially 
made when computing the log-intensity variance with photo- 
weaeis, is crucial in extinction measurement (which depends 
on accurate determination of the average detector voltage). 
The immediate concern was the possibility of a functional 
relationship between average detector voltage (V) and the 
exponential of the average log detector voltage (exp (In (V) ) ) 
dependent on the average value (l1n(V)) and the standard 
deviation (co) of the lognormal distribution. 

A computer program was developed which analyzed the 
variation oN the ratio of the two quantities with increments 
of either In(V) oro. The method was to assume an l1n(V) 
value for a Gaussian distribution and compute the average 
voltage. The probability density function of ln(V)was 


assumed, therefore, to be of the form: 


(x - xX) 
Y = C exp[- ] 
in 20 


where: 
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x = In(V) 
x = ln(v) 


Let y(V) be the probability density function of V: 


b exp (bd) 

f y1,9* = i y(V) av 
a exp (a) 
fe 3 Shatin we o 

: a dv Yin 
sine Yao = aye CV 

fee) dn (3) Yin fee) 

if Yint* = if T_ av = f y (V) av 

—< 1n( 0) 0 


Since the number of total counts accumulated over a range 
of detector voltage while sampling will be identical whether 


the display distribution is linear or logarithmic: 


Yin 
y (V) av = Yind% — inn av 
. Yin (x-x) 7 GC (1n (Vv) -In(v) * 
-» y(V) = 7 al al ayy, SS SS 
fe) FAQ} 
Applying statistics: 
Tye) GV. 


VY = 
fy (v) av 
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For a normal distribution a dispersion index of +30 gives 

a confidence limit of 99.73% [Ref. 7, et.al.) and these 
limits were imposed in a program applying Weddle's rule 
[Ref. 8] to accomplish the integration. The ratio of V 

to exp(ln(V)) was found to be independent of In(V) but an 
increasing function of o. With points generated by the pro- 
gram (for o values between 0.1 and 1.2 at increments of 0.02) 
a fourth-degree polynomial was found (using least-squares 
fitting) to generate average detector voltage values from 
average log-intensity. The polynomial coefficients were 
entered on the magnetic card carrying program constants 

and used to find the average detector voltage after deter- 
mining In(V) (from the coefficients of the least-squares 
polynomial fitting data to a lognormal distribution). An 
easier method was subsequently employed which can be 


justified by proceeding further with the preceding derivation. 


y= Sepnay 
but 

V = exp(ln(V)) = exp(x)- 
and 


y(v) dv = Vian ax 
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—, 2 
fC exp (x) exp [-2-%)_ 1} ax 


7 = He xP KX) Yin?x = 26 
= (Spo ee (x-%) 7 
ileal fc exp [-——*—] dx 
20 


After combining exponentials and completing the square in 


the numberator, and defining x' = x - ae 
= ae (x-x') 7 
C exp (x+=>-) fexp [-——5——] dx 
= 20 
=a 
a sexp [3d] ax 


2G 


The integrands in numerator and denominator are Gaussian 
distributions with the same standard deviations and, 


therefore, the same areas. 
a _ 2 ZA 
~ Ve= exp (x + ae eo (li) exp (5) 


This derived result can be found in meticulous statistical 
literature [Ref. 9]. Since ln(V) ando were already avail- 
able in the program the necessary calculation was easily 
included and required significantly less memory than the 
fourth-degree polynomial conversion. Incidentally, it was 
found that the polynomial method was accurate within 23% 
for standard deviation values between 0.1 and 0.9! 

With an assumed extinction constant the calculation of 
an extinction coefficient is easily made, i1.e.: 

Ln (>) 
VR 


R 
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where: 


uw = extinction coefficient 


C = extinction constant (VoR~ ) 


<l 
ll 


average detector voltage 


R = range 


This computation was added to the program with the extinction 
constant becoming one of the operator input parameters. 

The next program modification was the addition of a 
section which could process linear data to find the average 
detector voltage and extinction coefficient. The main 
purpose was to correlate results with those from the log- 
arithmic section but this was not done to any great extent. 
The average channel number is determined by dividing the 
summation of channel number times channel counts by the 
Summation of counts. The average channel number is then 
modified in a linear scaling process to give the average 
voltage. The extinction coefficient was computed in exactly 
the same way as in the logarithmic section. 

After the linear portion was added the program almost 
filled available memory and was found to execute more slowly 
than before — especially plotting data. The program showed 
minimal utilization of subroutines and, consequently, the 
program wasted memory. The subroutines used were scattered 


throughout the body of the program which meant wasted time 
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spent searching after branching statements. The computer 
still took data input only from the PIP-400 through the 
interface. A Nuclear Data 128 channel analyzer which inter- 
faced with a teletype was available as a back-up but the 
program had to be modified to accept manual input anda 
different channel capacity. With all the discrepancies 

a major reorganization was needed. | 

The HP981LO0A calculator branches to subroutines by going 
to the program beginning and searching through until finding 
the called label. This requires an inordinate amount of 
time to complete "do-loops" and subroutines. This problem 
was solved by making processes which were reiterative or 
used by more than one section of the program into subroutines 
and placing them in front of the main body. "Flag" state- 
ments were used to allow automatic or manual data input from 
the PIP-400 or Nuclear Data analyzer respectively. 

The original program used a large amount of memory in a 
ealeculation to filter out bias introduced in the scintillo- 
meter. When the single-detector system for a chopped laser 
source was first put together the demodulator would output 
a few tenths of a volt with no input present thus biasing 
the sample distribution. This offset had to be taken into 
account when processing data because of its effect when 
using the logarithmic scale. In subsequent demodulators 
the offset problem was solved and the calculation became 


extraneous. When the process was removed from the program 
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the only discernable difference was an increased speed in 
plotting data. 

The net result of all the program modifications pre- 
viously discussed was a drop in the average processing 
time of an experimental distribution from over ten minutes 


to less than six minutes. 


B. METEOROLOGY PROGRAMS 

Just prior to the most recent R/V Acania experiments 
(May, 1976) Professor Schacher (associated with the atmos- 
pheric turbulence group) requested a computer program to fit 
a Gaussian curve to data manually input from the Nuclear 
Data pulse height analyzer. No explanation of the relevant 
theory was given or required to distill the necessary 
processing from the existing scintillation/extinction program. 

Upon analyzing the processed data it was found that the 
Gaussian fit was a poor one. Prof. Schacher suggested the 
resultant data distributions appeared more Lorentzian and 
requested another program fitting data to that distribution. 


The Lorentzian distribution is of the form: 


Yor a 

ré -- ine) <4 ee oe ne 
- a = dt x2 = O58 x + eels a is 
[7 ¥ A A A 


This can again be fitted by least-squares method to a 


quadratic and the coefficients evaluated to get the average 
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voltage (after scaling x, the average channel number) and [ 
(after scaling). The weighting scheme differed from that 
used for the Gaussian. The uncertainty is (again assuming 


the number of counts in each channel is Poisson-distributed) : 


iL -2 -d -vy -3/2 
oe) te 8 dy = Sh = SY = -y 9 
i Y 
Again the data was weighted inversely with its uncertainty 
(multiplied by the number of counts in the channel to the 


3/2 power). The least-squares fit attempts to minimize: 


rw (Jae - bx. - c) 
aL Si a 
i 
where 
_ 3/2 
Op YG 
N/a = number of counts in channel 1 
x = channel number i 
ee = 
A 
_ -2x 
SS aR 
z _ x? +7? 
A 


Since all the machinery of the least-squares fitting of 


a quadratic applying Cramer's rule already existed, it was 
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a relatively simple matter to modify the existing Gaussian 
program. The comparative results definitely point to a 


Lorentzian distribution, which was not predicted by theory. 
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Sad, A SIMPLIFIED FLOW CHART 
ENTER MONTH, DAY, YEAR OF THE NEW PROGRAM 







NES 





CALIBRATION 
POINT PRINTOUT? 





CALIBRATION POINTS 
PRINTED OUT 


CONSTRUCT SCALE FACTOR 


ENTER CALIBRATION CHANNELS, 
CALIBRATION VOLTAGES 


ENTER EXTINCTION CONSTANT, 
WAVELENGTH, RANGE, PAR GAIN, 


MAXIMUM COUNT, TIME 
















LINEAR OR _INEAR 
LOGARITHMIC VOLTAGE 
DISTRIBUTION? 
IC ENTER DATA, CONSTRUCT SUMMATIONS 
Boos Mall FOR AVERAGE VOLTAGE 

PLOT AND LABEL AXES MWERCRo aon) OYE 

ENTER DATA AND PLOT. CON- 
CONSTRUCT AND PRINT EXT MU 
STRUCT SUMMATIONS FOR LEAST- 
ere eT. (EXTINCTION COEFFICIENT) 
CONSTRUCT COEFFICIENTS OF [END > 


LEAST-SQUARES QUADRATIC 


CONSTRUCT LOGARITHMIC 
VOLTAGE VARIANCE 


PLOTTER PRINTS DATE, TIME, 
BD (WAVELENGTH), RNGE 
RANGE), $ SQ (VARIANCE), 
STANDARD DEVIATION) 


CONSTRUCT AND PLOTTER PRINT 
EXT MU (EXTINCTION COEFFICIENT) 
PLOTTER PRINT GAIN (PAR GAIN) 









> PLOTTER PLOT COMPUTED GAUSSIAN 


CURVE AND CALIBRATION MARKERS 





CONSTRUCT AND PLOTTER PRINT 
CNS@ (C, SQUARED), CN ic) | END 


CONSTRUCT AND PLOTTER eng 
AVE V (AVERAGE VOLTAGE 
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001034 
000000 
000000 
000000 
006291 
000000 
000000 
000000 
000000 
006638 
000000 
000000 
000000 
000000 
000000 
000000 


*000000 
*000000 
#000000 
*000000 
*008860 
*000000 
*000000 
*000000 
#000000 
#003506 
*000000 
*000000 
*000000 
*000000 
*000000 
*000000 


*000001 
*000000 
*000000 
*000000 
*005971 
*000000 
*000000 
*Q00000 
*000000 
*Q00000 
*000000 
*000000 
*Q00900 
*Q000000. 
*000000 
*Q000000 


*000000 
*000000 
*000000 
*000000 
*000000 
*Q000000 
*QQ0001 
*000000 
*000000 
*000000 
*000000 
*006799 
*Q00000 
*000000 
*000000 
*000000 


*000000 
*000000 
*000000 
*000000 
*Q000000 
*000000 
*009367, 
*Q000000 
*000000 
*000000 
*000000 
*OO2LI922 
*000000 
*Q00000 
*000000 
*000000 


*000001 
*000000 
*000000 
*000000 
*000000 
*000000 
*016270 
*000000 
*000000 
*000000 
*000000 
*000000 
#000000 
*OOOL 2 i 
*000000 
*Q00000 


*000000 
*000000 
*000000 
*000000 
*000000 
*Q00000 
*000002 
*000000 
#000000 
*000000 
*000000 
*000000 
*000000 
*007691 
*000000 
*000000 


*000000 
*000000 
*Q000000 
*Q000000 | 
*000000 
*000001 
*000000 
*000000 
*000007 
*Q000000 
*000000 
*000000 
*Q00000 
*003369 
#000000 
*000000 


An example of the teletype output of the voltage calibration — 


channel numbers used in scaling 


(reproeduced from 1 March 1976). 


Notice the extraneous counts collected in channel l. 
points would be entered as follows: 


Channel # 


Voltage 
02275 


Jee 
LT. 0© 
2. 00 
4.00 


627097 
S203 
ZS 
2g 1.50 
1 EORS 0 


Change in # 


Lovo 
oa 7 2 
Be yw 
JE LEE 


The 





PICTURE OF EXPERIMENTAL DISTRIBUTION AS DISPLAYED ON THE 
PIP-400 PULSE HEIGHT ANALYZER 


49 





WWuoodd AIO AHL WONT LNdLNO YaLLOTd JO WIdWVWxXd NW 


WAZ AS) AS!) Ab! Az) A) (3S eS eh ee 





a> BAZ Se NS ‘eel —- LSS. = LSS 
BAe a = HEWSS iSama= 2GSaysenyps BWNoIsS 
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WWaS0dd MAN HHL WOU LBAdLNO YWALLOTd AO dIAMNVXad NW 


a.) -.1 ad | — Oe =| — i - L e 1 - 1 - 


°*L— 


* ik 





A AA fh 

“Aa | 

‘E—- Mes’! = MW LaLa = 
mai’: = ND> ‘hi- SSE! 
hes = = Ssa'2 

‘SEES = SNe rah val Pal = be ra 
"SSS1 BWIL “a2. 7°17 


Mun ioy 


+ Nipifa 
N AAH 
NSN> 
ns & 
daw) 
ae EU) | 


Dui 





SET FLAG FOR 


MONTH? 


CAL PRINIQUT 4, 
117012073 DAY? 
Bie 7 30. 
7040.103 YEAR? 
10768.104 Toe 
10883.134 FLAG - CAL PTS 
10939.164 PSE CAL? 
8756.194 61.0 
2194.195 2ND CAL? 
HNEPER 1ST CAL Lal YE, 
ENTER 2ND CAL 3RD CAL? 
ENTER 38D CAL 17 Bya0 
ENTER Vt IN Y V1l-Y AND V3-X. 
AND V3 IN X 0.50 
ENTER. VOLTAGE 8.00 
OFFSET. EXTINCT CONST? 
ENTER. WAVELENGTH 15300.00 
IN MICRONS WAVELENGTH 
ENTER. BANGE (MICRONS ) ? 
IN METERS 0.6328 
ENTER. DALE RANGE (METERS) ? 
ENTER TIME 1625. 
SEL ANALYZER TO GAIN? 
PRINT OUT. AND 200. 
SET FLAG FOR MAX COUNT? 
DATA PRINT OUT 8000. 
ENTER MAX COUNT TIME? 
13a 12 


FLAG - LINEAR 
FLAG = LOAD DATA 





DALE Pole PDN? 
20776. END. 
TIME 
4245. 
WAVELENGTH MONTH? 
O2O325 4, 
RANGE DAY? 
900. 28; 
YEAR? 
183.092 3 oe 
291.093 DSIRE OP Ny, 
503.094 2ND CAL? 
569.095 3RD CAL? 
645.096 Vi-Y AND V3-X% 
640.097 MAX COUNT? 
579.098 80000.00 
758.9099 TIME? 
714.100 HL pib Isic O20 
fA Ol LOAD DATA 
849.102 END 
etc 


Reproductions of computer printer tapes made during data 
processing: original program with cal and data printout 
(left), present program (upper rt.), meteorology program 
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Lid NES Smes 

2a2n (| 7eeSe Sa - A aA 
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APPENDIX A 


The following lists the connections made in mating the 
output plug of the PIP-400 to the input plug of the interface 
for the HP9Y81lOA: 


CALCULATOR 
LEAD FUNCTION PIP LEAD # LEAD # COLOR 
DATA 5 
roe = bre i 1 6 red 
2 2 7 white/black 
TA 3 8 white/green 
4 a oS 4 9 blue 
LO. = i. bit 5 10 white/blue 
20" 6 aval green 
-" 7 12 white/gray 
oa 8 ales yellow 
9-12 N/C 
ADDRESS 2 
LOe — le bee 13 30 white 
Dm 14 Sul white 
1 1esy, UGS: N/C 
Oe) =) 1 bate 17 34 white 
2 18 25 : 
4 19 36 C 
0 3 = 20 37 " 
10° = 1@it 21 38 white 
a. 22 39 : 
A PAS 40 : 
Sa 24 41 i 
DATA 3 
TO2"=—= lebit 25 14 white/violet 
—_ 26 ics white/black (2) 
4 @& a7 16 white/yellow/black 
5 oa. 28 ae. white/blue/black 
ite o= iebase 29 18 white/red/black 
2 (tm 30 1, white/orange/black 
4 "* Sj3e 20 black 
l S es 32 zie white/gray/black 
Loy = lL ore 33 aie white/yellow 
2. = 34 Ze white/brown/black 
4 " 35 24 white/orange 
0 oa. 36 25 white/green/black 
Oy = deepak 3 26 white/violet/black 
2° 38 27 white/red 
Aes ca 39 28 white/brown 
Bu bie 40 29 white/black 
Print Command 41 inverted 49 white 
Hold Command 42 ih 47 " 
Ground 43 chassis and power supply white 
44-48 N/C 
+ 6 volts 49 hex inverter white 
50 =. N/C 


N/C designates leads which are not connected 
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APPENDIX B 


The following is a listing of the data storage registers 
utilized by the program: 


REGISTER 


USE 


computed average logarithm of detector voltage 


computed summations for least-squares fit; 001, 
002 and 003 are used subsequently to store the 
coefficients of the least-squares polynomial 


9919 


determinants (formed from summations) which give 
coefficients of the least-squares polynomial 


computed standard deviation of least-squares 
Gaussian fitted to data 

first calibration channel number 

second I tt tt 

third te DY et 

computed variance of least-squares Gaussian, 
later used again in Cy Calculation 

time 

wavelength in meters 

range in meters 

wavelength in microns 

used in C,, calculation 


month s 
day 

year 

5/max count* 
gain 


extinction constant 

computed average detector voltage 

computed extinction constant 

data point 

data point weighting 

channel number-to-voltage scaling constant 
do-loop accumulators for the linear calculation 
of average detector voltage 

first calibration voltage 

chad is 4 

varying pen coordinate value for axis labelling 


incremental pen coordinace values for axis labelling 


label values in axis labelling 

236 (range of x-axis) 

11.8 (range of y-axis) 

alpha scale factor (0.1, FMT 11) 

computer-scaled x-coordinate 
y-coordinate 
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REGISTER 


Eis. * 
ECGs 
tO * 
0S = 


aé&b 


-0.9 OG wy ) 
10.9 Fane (followed by FMT 13) 
-18 (xX 7%) 
218 (xM1n) (followed by FMT 12) 


used many times through out the program 


* marks those registers whose values are input by the operator 
** marks those registers with values stored on a magnetic card 


The following is a list of label titles in the order 
they appear in the program: 


LABEL 
A* (UP, DN) 


Bemtel , 2) 


ees , 4) 


B2(A,B Ge, 5) 


FUNCTION 
takes manual or interface input; puts data in 
060, channel in a, 5/max count in b 
puts channel number in y and upper data channel 
arate xx 
puts channel number in y and upper channel 
Lane) x 
plots data points and constructs summations 
necessary for least-squares fitting 


plots the computed Gaussian and calibration points 


plots marker at calibration point 
constructs a special floating point form of 
a number to be plotted 
labels plotted x & y axes 
computes the extinction coefficient 
prints operator instructions: 

FLAG — LOAD DATA; PIP ON PRINT? 
lists calibration points 
preliminary program set-up (month, day, year, 
calibration points, extinction constant) 
completes program set-up (wavelength, range, 
gain, max count, time); branches to label N 
for linear; plots and labels x & y axes; plots 
data; calculates desired constants; plots con- 
stants, computed Gaussian, and calibration 
points; returns to label M 
linear portion of program; computes average 
detector voltage and extinction constant, 
and prints them, returns to label M 


* marks labels which are subroutines 
() hold subroutines and labels branched to 
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APPENDIX C 


The following is a complete listing of the program as 
printed out by the calculator; 


0000 GTO a4 0050 LBL si 0100 aE 23 
0001 LBL au 0051 DN 2 0101 DN 25 
0002 L i2Z 0052 LBL jo) 0102 CNT 47 
0003 LBL oye 0053 Be 27 Oro 3 S7R 77 
004 A 62 0054 FMT 42 0104 LBL Sys 
0005 IFG 43 0055 al O1 0105 B 66 
0006 SFL 54 0056 8 10 0106 a i 
0007 GTO 44 0057 0 00 0107 UP 27 
0008 LBL om 0058 1 O1 0108 IFG 43 
0009 UP Zi 0059 xXTO INS CEOS Wis oy 54 
0010 FMT 42 0060 a es: 0110 GTO 44 
0011 il Ol 0061 0 rs Offi sb alk 
0012 8 iO 0062 0 13 OZ 1 O01 
0013 3 03 0063 ne PLY OE s.> CNT 47 
0014 FMT 42 0064 CLX ad 0114 i O1 
0015 3 03 0065 STP 41 OE > 2 jib 
0016 3 03 0066 xXTO 23 0116 5 05 
0017 F oy 0067 6 06 Oi? Gro 44 
0018 We 27 0068 0 00 0118 LBL Syl 
Ong CLx a7 0069 FMT 42 Out? 2 02 
0020 xX=Y 50 0070 ib Ol 0120 LBL oi 
0021 GTO 44 0071 8 10 O21 al O1 
0022 LBL onl 0072 ) ONS: 0122 CNT 47 
0023 A 62 0073 LBL Sel 0123 1 Ol 
0024 CNT 47 0074 DN ZS 0124 2 02 
0025 EEX 26 0075 XFR 67 H125 3 03 
0026 CHS 22 0076 5. 05 0126 LBL oul 
0027 3 03 0077 4 04 O27 2 02 
C028 Xx 36 0078 xXTO JS) Ol2s S/R 77 
0029 DN 25 0079 b 14 0129 LBL ol 
0030 Ue Ae 0080 1/X IL 7 0130 C 61 
0031 INT 64 0081 DE 27 OS a 3 
0032 - 34 0082 al Ol 0132 UP 27 
0033 xXTO 23 0083 0 00 iss iFG 43 
0034 6 06 0084 DIV ee. 0134 SFL 54 
O35 0 00 0085 XFR 67 OLS 5) 135ue, 44 
0036 EEX 26 0086 6 06 0136 LBL spk 
0037 | 03 0087 0 00 Ons 7 3 O3 
0038 X 36 0088 xX<Y a2 0138 1 O1 
0039 Z 02 0089 CNT 47 OL39 and 
0040 0 00 0090 S/R 77 0140 3 dal 
0041 0 00 0091 CNT 47 0141 GTO 44 
0042 xX<¥Y D2 0092 CNT 47 0142 LBL ol 
yo43 CNT 47 ONO eis UP Psi) 0143 é. 04 
0044 - 34 0094 a i) 0144 LBL = 
0045 CNT 47 0095 Be 27 0145 3 03 
0046 CNT 47 0096 EEX 26 0146 a O1 
0047 YTO 40 UO 7 3 Oe: 0147 2 02 
0048 0 eS 0098 DIV 35 0148 8 10 
0049 GTO $4 0099 DN 25 0149 LBL opt 
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0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0208 
0209 
0210 
OZE 
OZ 
GZS 
0214 
O25 
0216 
0217 
0218 
OZ1y 
0220 
0221 
O22 
0223 
0224 
0225 
0226 
OZ 27 
0226 
0229 
0230 
0231 
i232 
0233 
0234 
0235 
0236 
O237 
0238 
0239 
0240 
0241 
0242 
0243 
0244 
0245 
0246 
0247 
C248 
0249 


22) 


0250 
0251 
GZa2 
0253 
0254 
QiZ55 
0256 
O25 7 
0) 5 ave 
0259 
0260 
0261 
0262 
0263 
0264 
0265 
0266 
0267 
0268 
0269 
0270 
O271 
0272 
27.3 
0274 
0275 
0276 
0277 
0278 
0279 
0280 
OZ ok 
Oa 
0283 
0284 
OZES 
0286 
02 ay 
0288 
0289 
0290 
OZ91 
0292 
O293 
0294 
UZ 35 
0296 
0297 
0298 


(0299 
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0500 
OSOe 
0502 
0503 
0504 
0505 
0605 
0507 
0508 
0509 
0510 
O5i1 
Gad 
0513 
0514 
OSS 
O26 
0517 
0518 
0519 
0520 
0521 
OS22 
0S25 
0524 
0525 
0526 


0527 


0528 
Oa20 
0530 
Om sal 
0532 
OSS 
0534 
OBS. 
0536 
0537 
OFS) S32) 
esyey, 
0540 
0541 
0542 
0543 
0544 
0545 
0546 


0547 


0548 
0549 


Gul 
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0800 


0802 
0803 


0807 
0808 
0809 
0810 
O8ll 
0812 
O) 5 ls, 
0814 
0815 
0816 
0817 
0818 
0819 
0820 
0821 
0822 
0823 
0824 
0825 
0826 
0827 
0828 
0829 
0830 
O's I 
0832 
OSS 
0834 
O16 35 
0836 
Gus 7) 
0838 
0839 
0840 
0841 
0842 
0843 
0844 
0845 
0846 
0847 
0848 
0849 
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64 





65 





66 








67 





re io ~ ay 


p UX 
OZUNNHHARAWA.- 
<= © 


ze) 


68 
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1750 
IEFyoms 
Ey SZ 
33 
1754 
7S 
E76 
IES) 7 
JEY LS): 
1759 
1760 
OFS I. 
IAS 
1763 
1764 
ley sy) 
IL FASS) 
1767 
iG 
ive9 
1770 
77 e 
ec 
Da 3 
1774 
dF ey Re. 
1776 
er Tel 
IE TEs} 
ILE TPS, 
1780 
yok 
to 
i733 
1784 
1785 
1786 
v7.37 
1788 
Wipes, 
Lee 
eo 
IL 4 
eS) 
1794 
IL yes 
1796 
oe 
Igoe 
IoD 





Li) 8)0. 
1851 
TE52Z 
1355 
1854 
S55 
1856 
ILS) 7s 
1858 
i359 
1860 
1861 
1862 
1863 
1864 
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